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I extend the well-known photonuclear sum rule that relates the strength of the photoexcitation of 
the giant dipole resonance in a nucleus to the number of elementary scatterers-nucleons to the case 
of virtual photons. The new sum rule relates the size of the magnetic polarizability of a nucleus 
to the slope of the transverse virtual photoabsorption cross section integrated over the energy in 
the nuclear range. I check this sum rule for the deuteron where necessary data is available, discuss 
possible applications and connection with other sum rules postulated in the literature. 
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Scattering of light off a composite object has long been 
used to study its structure. At low frequencies, elec¬ 
tromagnetic waves scatter without absorption and solely 
probe its mass and electric charge, the classical Thomson 
result. With the photon energy raising above the absorp¬ 
tion threshold internal structure is revealed. Kramers 
and Kronig related the photoabsorption spectrum of a 
material to its index of refraction by means of a disper¬ 
sion relation mm based on the probability conservation 
and causality. Dispersion relations and sum rules have 
been among the main tools for studying the electromag¬ 
netic interactions in atomic, nuclear and hadronic physics 
domains. These domains roughly correspond to keV, 
MeV and GeV photon energies, respectively, and this 
scale hierarchy indicates that dynamics in each domain 
can be clearly identified. Thomas-Reiche-Kuhn sum rule 
equated the sum of oscillator strengths in an atom to the 
number of electrons [3HS]. For nuclei, Levinger-Bethe [5] 
and Gell-Mann, Goldberger and Thirring [7j related the 
integrated photoabsorption cross section to the number 
of elementary scatterers, protons and neutrons in a nu¬ 
cleus. For GeV energy photons that resolve the nucleon 
structure, Gorchtein, Hobbs, Londergan and Szczepaniak 
[5] observed that the integrated strength of the nucleon 
resonances may be explained by counting the constituent 
quarks. These sum rules are an economic, albeit approx¬ 
imate way to express duality, the transcendence of higher 
energy degrees of freedom in the low-energy phenomena 
[9]. In this letter I extend the Thomas-Reiche-Kuhn- 
Levinger-Bethe sum rule to the case of virtual photons, 
obtain a sum rule for the nuclear magnetic polarizability, 
and discuss further applications. 

The spin-averaged, forward Compton tensor is ex¬ 
pressed in terms of two scalar amplitudes Tl i2 (f, Q 2 ), 
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with the invariants defined in terms of the nucleus and 


photon four-momenta p, q as v = (p ■ q)/Mp , Q 2 = 
—q^q^ = —<? 2 > 0, and p 2 = Mj., with M T the target nu¬ 
cleus mass. In this letter I concentrate on the transverse 
amplitude T\ . Its imaginary parts is related to the unpo¬ 
larized structure function Fj as ImTi = (7ra ern /Mr)Fi, 
with a em ss 1/137 the fine structure constant. T) satis¬ 
fies a once subtracted dispersion relation (DR), 


ReTiiy, Q 2 ) = Ti(0, Q 2 ) + 


v-emV 2 f dv' 2 F\(y', Q 2 ) 
Mp J v' 2 (v ' 2 — v 2 ) 
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where the integral is understood in terms of its princi¬ 
pal value. I remove the pole contribution that is due to 
an absorption of a virtual photon by an on-shell ground 
state (this separation is well-defined, see, e.g., discus¬ 
sion in HD!)- Upon this removal, the subtraction con¬ 
stant T™ p (0, Q 2 ) is defined in terms of the nuclear charge 
form factor Fc normalized to unity at Q 2 = 0, and the 
nuclear magnetic polarizability P^ cl {Q 2 ) generalized to 
finite Q 2 , 
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with Z(N) the number of protons (neutrons) in the nu¬ 
cleus, a em « 1/137 the fine structure constant, M « 
M p ss M n the nucleon mass, such that Mp ~ (Z + N)M. 

Real photoabsorption on lead, shown in Fig. [l] il¬ 
lustrates several general features common to all nuclei: 
i) the strength of nuclear excitations is concentrated in 
the region between the breakup threshold UminiQ 2 ) = 
B + Q 2 with B the nucleon removal threshold 
for the nucleus, and v max (Q 2 ) « B + Q 2 /(2M) + 30 
MeV; ii) nuclear cross sections stay small above that en¬ 
ergy and below the threshold for the nucleon breakup 
VniQ 2 ) = Q 2 /(2M) + to t + m 2 /(2M), with m, the pion 
mass; iii) above this threshold, an incoherent absorption 
by Z protons and N neutrons that make up a nucleus is a 
good overall representation of the cross section (modulo 
nuclear effects). I exploit the observed gap between v max 
and by evaluating the DR for X) at an intermediate 
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FIG. 1: (Color online) Total photoabsorption cross section on 
lead in /ibarn as function of energy. Data in the nuclear range 
(blue crosses) extend up to v max ta 30 MeV are from Ref. 
nu. Data above the pion production threshold (red open 
circles) are from [12H15) . The vertical dashed lines display 
Vmax,i / oa, and , see text for further details. 


energy (Q 2 ) ss B + Q 2 /(2M) + 70 MeV, impose the 
hierarchy of scales, v 2 nax <C <C z/ 2 and take respective 
limits, 

V-max 

ReT r(^oo,Q 2 ) = T? p (0,Q 2 ) -F^Q 2 ) 

IS min 
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For compactness, I suppressed the Q 2 -dependence of 
the integration limits. The integral between v max and 
is understood in the sense of its principal value. Next, 
the scale hierarchy is used to calculate ReT" p (i/ oc , Q 2 ): 
the scale Vqc was chosen such that the bulk of nuclear 
excitations lies significantly below it. Then, photons will 
scatter off essentially unbound nucleons; the energy is 
significantly lower than the pion production threshold, 
so the nucleon structure is not resolved at that energy, 
and it is legitimate to approximate its value by a low- 
energy expansion up to order z/^, 

Re7T>oo,Q 2 ) = -Z^FtfiQ 2 ) - N^F^iQ 2 ) 

+ZQ 2 p p M (Q 2 ) + NQ 2 (3^(Q 2 ) 

+ W4 J°° dv q2) + NF n {Uj q2)] 5 (5) 

where F^'" 1 denotes the proton (neutron) Dirac form fac¬ 
tor, and /3^"^(Q 2 ) stand for the proton (neutron) mag¬ 
netic polarizability, respectively, extended to finite Q 2 . 
A subtracted dispersion relation analogous to that of Eq. 


(|2j) is imposed on the single nucleon amplitudes, with 
pj >,n free nucleon structure functions. Now, Eqs. m 
can be combined together, and the coefficients at differ¬ 
ent powers of equated. If nuclear and hadronic scales 
are indeed well-separated, above v m ax{Q 2 ) nucleons are 
unbound, and the coefficient at should vanish, 
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Turning to the terms independent of z^, and setting 
Q 2 = 0 Levinger and Bethe [5] obtained, 

r^ma.n: f l.. 

ZN = 2 -F\(zz,0), (7) 

J Vmin V 

i. e. integrated strength of nuclear excitations is fixed by 
the number of nucleons within the nucleus. Levinger- 
Bethe sum rule of Eq. Q is obeyed for a wide range of 
nuclei, typically better than 10% [IS]- As an example, 
the parametrization of the deuteron photodesintegration 
cross section in Ref. m leads to the value of the right 
hand side 1.007, in excellent agreement with the sum 
rule, NZ = 1. Deviations due to non-vanishing of the 
principal value integral and effects of nuclear binding and 
shadowing in Eq. © were estimated, e.g., in Refs. HID- 
I now consider the first derivative with respect to 
Q 2 at the origin. Using the charge radius defined as 
R 2 Jh = — 6F^(0), the sum rule for the nuclear magnetic 
polarizability is obtained, 
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where I neglected effects of nuclear and nucleon recoil 
that enter the Q 2 -dependence of the integration limits 
(above taken at Q 2 =0), effects of nucleon polarizabili¬ 
ties and nucleon charge radii. 

This sum rule is useful since for most nuclei the 
magnetic polarizability is not known, unlike the sum 
a nuci pnuci t kat j g g X ed by Baldin sum rule T8], 
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and can be directly extracted from the experimental data. 

To my knowledge, deuteron is the only nucleus for 
which theoretical predictions of f3^ cl exist, calculated 
in EFT [19] and potential model (20] approaches, sum¬ 
marized as /3^ = 0.072(5) fm 3 . One can now check, how 
important the neglected terms are numerically. Using the 
value of the proton charge radius from recent [iH mea¬ 
surements eh usi, and the neutron charge radius along 
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with the nucleon magnetic polarizabilities from the PDG 
[28] gives ~ 1.6 x 1CU 3 fm 3 , two orders of magnitude be¬ 
low /3fj. The effect of the deuteron charge radius taken 
from [23. is of the similar order, ~ —1.5 x 1CU 3 fm 3 , 
also negligible. However, for heavy nuclei these two con¬ 
tributions can have very different size, e.g., for lead the 
two terms give ~ 0.08 fm 3 and ~ —0.5 fm 3 , respectively, 
which explains the choice of keeping the nuclear radius 
effect but neglecting the nucleonic contributions. The 
value of Pm for lead is unknown, but cue + Pm ~ 14.5 
fm 3 [H] gives a rough idea, even though it can be ex¬ 
pected that Pm < 0.1 cxe for that nucleus. 

Using a recently proposed detailed parametrization of 
deuteron breakup data m that covers Q 2 in the range 
[0.005 GeV 2 ; 3 GeV 2 ] and energy between the deuteron 
breakup threshold and well into the hadronic range, a 
numerical evaluation of the right hand side of Eq. ([8]) 
can be done. It leads to P d M = 0.096(15) fm 3 , close to the 
model-based expectation, P^ = 0.072(5) fm 3 . Note that 
even raising v max to 140 MeV would increase the integral 
by mere 1%, so the result is very robust. To enforce the 
agreement, one needs to modify the parametrization of 
Ref. [17] (Eq. (27) and Table II of that Ref.) via 


to 


2.15(35) x 10 4 GeV" 3 Q 2 
(1 + 52(8) GeV~ 2 Q 2 ) 2 


fFSi (o2) = 1-61(11) x 10 4 GeV" 3 Q 2 
lT W> (1 + 35(6) GeV _2 Q 2 ) 2 - 2 


( 10 ) 


( 11 ) 


The error in the numerator is fixed by that in the value 
of P'm , and the error (and a different power) in the numer¬ 
ator is obtained by a new fit to the quasi elastic data, as 
described in Ref. H7j. The two fit functions are shown in 
Fig. [2] With this exercise I demonstrate that the existing 
deuteron quasi elastic data are consistent with the pro¬ 
posed sum rule. The original parametrization Ref. m 
led to a 1.5cr-disagreement because the slope parameter 
was obtained by an extrapolation beyond the kinemati- 
cal range covered by the data without using the value of 
Pm as a constraint. 

Another sum rule involving the Q 2 -slope of the inte¬ 
grated structure functions was proposed by Bernabeu 
and Jarlskog [25]. They assumed that the longitudinal 
amplitude obtained as a linear combination of T) and X 2 
obeys an unsubtracted dispersion relation, and argued 
that the longitudinal structure function has to vanish 
identically at the real photon point independently of the 
energy to ensure gauge invariance, hence the integral be¬ 
comes convergent. In this way they arrived at a sum rule 
for the electric polarizability cxe alone, which is however 
incompatible with the Pm sum rule proposed here. I be¬ 
lieve that the reason for the disagreement lies in their 
use of unsubtracted dispersion relation. Since it is the 
Q 2 -slope that gives the sum rule, one in reality explicitly 



FIG. 2: (Color online) The comparison of the old fit without 
the use of the sum rule, f£ SI (blue dashed curve) and the 
new fit using the sum rule, fj- SI (red solid curve), with the 
uncertainty of each fit indicated by the band of the respective 
color. The sum rule is indicated by the star. The shaded 
band shows the kinematical range that is covered by the ex¬ 
isting D(e,e')pn data. The subview in the upper right corner 
magnifies the small values of Q 2 where the the slope of the 
new fit function is fixed to reproduce the value of Pm- Data 
points correspond to experimental data sets analyzed in Ref. 
[17] (Refs. [35-42] of that article). 


departs from the real photon point; then, the argument 
of vanishing of the longitudinal structure function at in¬ 
finity is no longer valid, and one is left with a divergent 
integral, so that the limit Q 2 —>■ 0 does not exist. 

The parametrization of deuteron quasi elastic data was 
used in Ref. m to estimate the two-photon exchange 
correction to the 2 P — 2 S Lamb shift in muonic deu¬ 
terium atom. A modification of the data parametrization 
proposed above based on the new sum rule will lead to 
a different prediction for that correction. Moreover, the 
photonuclear sum rule discussed above can further be ex¬ 
tended beyond its value and slope at Q 2 = 0 (TRKLB 
and the P™' d sum rule, respectively) to predict the full 
<5 2 -dependence of the subtraction function via 

Tr P (0,Q 2 )-T” p (0,0) (12) 

iw*(Q 2 ) 

= 2 W / V [F,M 2 )-FiM)], 
*w(Q 2 ) 


which contributes to the shift of the 2 S state through 


OO 

AE%s b = Aaemplsp))J d Qll(Tl) 


TH0,Q 2 )-Tr(0,0) 

Q 2 ’ 


(13) 


with 7 i(a:) = (1 — 2x)^/l + x + 2a: 3 / 2 , 77 = Q 2 /(4m 2 ), 
mi the lepton mass, and c^ 2 s (0) = (Za em m r ) 3 /irn 3 the 
squared atomic wave function at origin with the reduced 
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A E\ s 

This work 

Ref. Q71 

Refs. ETH50] 

A 

-2.294(740) 

-2.357(740) 

- 

A E$g bt 

0.505(35)(40) 

0.763(40) 


A 

-1.945(740) * 

-1.750(740) * 

-1.709(15) 


TABLE I: TPE contributions to the shift of the 2 S state in 
muonic deuterium in units of meV. The inelastic contribution 
is a sum of “PWIA”, “FSI”, “_L” and “hadr” contributions 
listed in Table I of m- The numbers in the first and second 
column in this row correspond to the use of fx SI and /t S 7 , 
respectively. Subtraction contribution is calculated with the 
sum rule in this work, while the number in the second col¬ 
umn is a sum of “Th.” and “/?” terms in Table I of [17]. The 
total contribution is obtained by adding the upper two num¬ 
bers with elastic term obtained in El, and the star indicates 
the inclusion of the internal Coulomb correction of 0.261 meV 
[27] . Total contribution summarizing potential models calcu¬ 
lations m is listed in the rightmost column. 

mass m r = M T mi/(M T + mi). The value of 7\(0,0) is 
subtracted to account for its inclusion in the lowest order 
atomic calculation. A similar approach based on the fi¬ 
nite energy sum rule obtained upon removing the Regge- 
behaved part of the hadronic photoabsorption, was ap¬ 
plied to the muonic hydrogen Lamb shift [28]. 

Numerical evaluation leads to a new prediction (in¬ 
cluding the inelastic Coulomb correction, as pointed out 
in Ref. [27]) A = — 1.945(740)(40) meV, to be com¬ 
pared to A E™ = —1.750(740) meV [17] . The two values 
agree within the error that is dominated by the uncer¬ 
tainty due to the low-Q 2 behavior of the quasi elastic 
cross sections. The systematical uncertainty in the sec¬ 
ond bracket is due to the use of the sum rule for the sub¬ 
traction term, and was estimated by varying the value 
of Vmax between 30 MeV above the quasi elastic peak, 
and the pion production threshold. An additional 0.01 
meV uncertainty due to was added in quadrature. 
It amounts in « 8% uncertainty and can be compared 
to 1% in the sum rule for The reason for the larger 
uncertainty is mostly in a steep rise with Q 2 of the QE 
peak that resides at higher energy than the threshold 
peak that completely dominates at Q 2 = 0. 

The large uncertainty of the DR result at present 
prevents one from talking of a disagreement between 
the new prediction and other models, nevertheless when 
new deuteron quasi elastic data at lower Q 2 will be¬ 
come available [31] the uncertainty may be sizably re¬ 
duced e- In that case the shift of —0.195 meV will 
result in a different value of the deuteron charge radius 
extracted from the iD Lamb shift measurement. Us¬ 
ing A E^g = 6.1103(3)(i?d/fm) 2 meV [27], the extracted 
value of Rd would be larger by SRd = 0.007 fm. It is 
smaller than the uncertainty of the radius extraction from 
scattering data R^ D = 2.128(11) fm but considerably 
larger than that using the isotope shift measurements 
[321 33] and muonic hydrogen Lamb shift EH ED , as 


well as the expected uncertainty of the muonic deuterium 
data. The method based on the new sum rule provides 
a different basis for estimating the subtraction function, 
as compared to the minimalist assumption used in Ref. 
El that the Q 2 -dependence of the deuteron magnetic 
polarizability resembles that of the charge form factor 
Pm(Q 2 ) ~ Pm^c(.Q 2 )- The sum rule-based calculation 
can be seen as a valuable systematic study of DR calcu¬ 
lations. A direct calculation of /?^(Q 2 ), e.g., in an EFT 
approach would help further assessing this systematics. 

The method proposed here can be used for calculating 
the subtraction function contribution to the Lamb shift 
in other light muonic atoms with the new experiments 
underway 1341 . For nuclei beyond deuteron a reliable es¬ 
timate of in potential models and in effective the¬ 

ories might be considerably more complicated. The pro¬ 
posed sum rule may serve a model-independent tool to 
extract from data, e.g. interpret measurements of 

Ml strength in heavy nuclei [351138] , 

Currently, models of a strongly bound composite Dark 
Matter (DM) [37 a have received much attention. Such 
DM particles would have electromagnetic polarizabilities 
and could interact with ordinary matter by means of the 
two-photon exchange [35]. At present, estimates of the 
nuclear part of the interaction have a modest ± order of 
magnitude accuracy [38]. For more quantitative calcula¬ 
tions based on dispersion relations the new sum rule will 
help constraining the subtraction function contribution. 

In summary, I proposed a new sum rule that general¬ 
izes the Levinger-Bethe sum rule to the case of virtual 
photons. Its slope at zero photon virtuality relates the 
nuclear magnetic polarizability to the slope of the trans¬ 
verse photoabsorption cross section integrated over the 
nuclear energy range. I showed that the quasielastic data 
on the deuteron are compatible with the sum rule, and 
applied its full version to the calculation of the Lamb 
shift in muonic deuterium. I discussed applications to 
light muonic atoms and direct DM detection. 
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